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Abstract 



OO ' We evaluate the Schechter-Valle (Black Box) theorem quantitatively by considering the 

most general Lorentz invariant Lagrangian consisting of point-like operators for neutrino- 
less double beta decay. It is well known that the Black Box operators induce Majorana 
<~i \ neutrino masses at four-loop level. This warrants the statement that an observation of 

Qh' neutrinoless double beta decay guarantees the Majorana nature of neutrinos. We calcu- 

late these radiatively generated masses and find that they are many orders of magnitude 
O ' smaller than the observed neutrino masses and splittings. Thus, some lepton number 

i-^ \ violating New Physics (which may at tree-level not be related to neutrino masses) may 

induce Black Box operators which can explain an observed rate of neutrinoless double 
P\| \ beta decay. Although these operators guarantee finite Majorana neutrino masses, the 

^ ' smallness of the Black Box contributions implies that other neutrino mass terms (Dirac 

or Majorana) must exist. If neutrino masses have a significant Majorana contribution then 
this will become the dominant part of the Black Box operator. However, neutrinos might 
\ also be predominantly Dirac particles, while other lepton number violating New Physics 

■ dominates neutrinoless double beta decay. Translating an observed rate of neutrinoless 

(<~^ \ double beta decay into neutrino masses would then be completely misleading. Although 

the principal statement of the Schechter-Valle theorem remains valid, we conclude that 
the Black Box diagram itself generates radiatively only mass terms which are many orders 
of magnitude too small to explain neutrino masses. Therefore, other operators must give 
the leading contributions to neutrino masses, which could be of Dirac or Majorana nature. 
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1 Introduction 



During the last decades neutrino flavor transitions were eventually found in oscillation exper- 
iments with atmospheric, solar, reactor, and accelerator neutrinos. These experiments have 
shown without doubt that neutrinos have small non-vanishing masses [1~3] requiring some 
extension of the Standard Model. Different possibilities exist to generate neutrino masses, 
but their special quantum numbers allow them to be Dirac or Majorana particles. This issue 
is directly related to the question whether lepton number is or is not a symmetry of Nature, 
because a Majorana mass term for the neutrino violates lepton number by two units. This 
question cannot be answered by standard neutrino oscillation experiments, unless one would 
be able to directly detect a conversion between neutrino and antineutrino (or vice versa). 
Therefore, other experiments have to be performed to determine the nature of neutrinos. 

Unfortunately, lepton number violating processes generically have very small amplitudes, 
as they are usually suppressed by tiny neutrino masses. Therefore, it is very difficult to 
observe these processes experimentally. At the moment, the most promising attempts to find 
lepton number violation are the experiments searching for neutrinoless double beta decay 
{Ou/3/3). Several such experiments have been performed (IGEX [4], Heidelberg-Moscow [5], 
CUORICINO [6], NEMO [7], and others), so far without unambiguous detection. New second 
generation experiments like GERDA [8], which has been built and is being readied for data 
taking, will soon improve the sensitivity and may observe We would like to emphasize 

that such a signal first of all establishes lepton number violation and that one should be a bit 
careful not to translate it immediately into a Majorana neutrino mass. We will therefore take 
a closer look at the well-known Schechter-Valle theorem [9], which guarantees that neutrinos 
are Majorana particles if neutrinoless double beta decay is observed. We will determine 
the numerical magnitude of these guaranteed Majorana mass contributions and discuss the 
consequences. 

The paper is organized as follows: We first review the classic version of the Black Box 
theorem, as well as its extension to general lepton number and lepton flavor violating processes 
in Sec. 2. In Sec. 3, we classify all point-like operators that can contribute to 0^/3(3 and 
calculate the mass corrections for two examples of operators, one of which leads to a vanishing 
correction while the other one generates a tiny but non-zero neutrino mass correction. We 
finally explain the connection of our results to chiral symmetry breaking in Sec. 4 before we 
conclude in Sec. 5. 

2 Schechter-Valle Theorem (Black Box Theorem) 

The Schechter-Valle theorem (Black Box theorem) [9] relates the effective neutrinoless double 
beta decay operator to a non-zero effective Majorana electron neutrino mass. The underlying 
argument can easily be explained using Fig. 1. Effectively, neutrinoless double beta decay 
may be seen as a "scattering" process (0 — )• dduuee) because, under the assumption that 
the weak interaction is described by a local gauge theory, crossing symmetry holds. One only 
has to additionally assume that [10] 

1. the u and d quarks and the electron are massive, and 

2. the standard left-handed interaction {T^^^cl + uZ^^dL) exists. 
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Figure 1: Contribution of the Black Box operator to the Majorana neutrino mass [9]. 



Then, it is possible to draw the diagram in Fig. 1, so that neutrinoless double beta decay 
induces a non-zero effective Majorana mass for the electron neutrino, no matter which is the 
underlying mechanism of the decay. The Black Box is nothing but an effective operator for 
neutrinoless double beta decay which arises from some underlying New Physics. The first 
assumption is necessary to ensure that two identical neutrinos are created. This can be seen 
in the following way [10]: We do not know anything about the chirality of the electrons and 
quarks produced by neutrinoless double beta decay. However, this assumption guarantees 
that we can make the particles running in the loops in Fig. 1 left-handed, by mass insertion 
if necessary. Thus the standard left-handed interaction from the second assumption produces 
the same type of neutrino at both vertices. Otherwise it would be possible that a neutrino 
and an antineutrino are created, which would give a Dirac mass term. 

Note, however, that the diagram in Fig. 1 is certainly not the only one that generates a 
non-zero effective Majorana mass for the electron neutrino. Other tree and loop diagrams 
exist and in addition the physical neutrino masses depend also on Dirac mass terms. Further- 
more, there may even be cancellations between different Majorana contributions which are 
induced by the Black Box diagram(s). This may appear as a fine-tuning, but the observed 
fermion mass patterns suggest that symmetries which explain these patterns may exist, and 
such symmetries could also lead to non-trivial cancellations. Taking into account this possi- 
bility of cancellations, Takasugi [10] and Nieves [11] improved the argument of Schechter and 
Valle [9], and showed that there cannot be a continuous or discrete symmetry protecting a 
vanishing Majorana mass to all orders in perturbation theory. We will follow the arguments 
of Takasugi [10] here. He assumed an unbroken discrete symmetry protecting the Majorana 
neutrino mass (the r/'s are global phase factors): 

T^eL rj^l^eL, ^eCL, QL ??ggL [Q = U,d), W^^ r]wW^^ . (1) 

To forbid the Majorana mass term, we need to have 

Vl^l, (2) 
and the invariance of the left-handed interaction requires 

vtVeVw = nLVdVw = 1 • (3) 
However, the existence of Oz^/3/3 (that is, the process di + ^ ul + ul + ei + cl) implies 

Vlr]frie='i-- (4) 

It is easy to see that Eqs. (2), (3), and (4) cannot be solved simultaneously. Thus, if 
the Majorana mass term is forbidden by an unbroken discrete symmetry, there will be no 
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Figure 2: Contribution of the entry of the Major ana neutrino mass matrix to the effective 
lepton number and lepton flavor violating vertex Fq,^. 




Figure 3: Contribution of the effective lepton number and lepton flavor violating vertex F^/s 
to the entry of the Majorana neutrino mass matrix. 



neutrinoless double beta decay. On the other hand, if neutrinoless double beta decay exists, 
there cannot be a symmetry protecting the Majorana mass and the term will be induced (the 
possibility of accidental cancellation to all orders in perturbation theory is not excluded, but 
appears to be very unlikely). 

The implications are quite strong, and this becomes particularly clear when we think 
of contributions to neutrinoless double beta decay from New Physics beyond the Standard 
Model. Besides the usual mass mechanism involving the internal exchange of a virtual Ma- 
jorana neutrino, many models of New Physics allow neutrinoless double beta decay also via 
other lepton number violating operators. Thus the theorem states that it is not possible to 
construct models with massless neutrinos where occurs. 

However, it cannot be overemphasized that the theorem so far is only a qualitative state- 
ment and does not say anything about the size of the induced Majorana neutrino mass. A 
quantitative determination of the mass generated for different realizations of the Black Box 
by the diagram in Fig. 1 is therefore a very interesting question. We will do such a calculation 
for different operators using an effective field theory approach independent of the underlying 
model in Sec. 3. The results will be interesting, as they clearly show the limitations of the 
Black Box theorem. 

Note that there exists also an extended Black Box theorem [12] which takes the three- 
generation Majorana neutrino mass matrix into account and includes the fact that i^e, t'^, 
and Vt mix. Moreover, the theorem has been extended to arbitrary lepton number and lepton 
flavor violating processes. It was also found that there exists a general set of one-to-one 
correspondence relations between the effective operators generating these processes and the 
elements of the neutrino mass matrix. In particular, AL = 2 processes described by ~^ la^p 
conserving baryon number have been discussed, which means that lepton number violation 
manifests itself via two external charged leptons in the final state. is a set of external 
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particles with B = L = and electrical charge Q = —2. Note that, for a lepton number 
violating process with the same flavor states la and different sets of external particles 
are possible. With a symmetry argument, similar to the one used before to prove the classic 
Black Box theorem, it has also been shown that the following relation between the effective 
lepton number and lepton flavor violating vertex Fq,^ and the entry of the Majorana 
neutrino mass matrix exists: 

M^^ = O^Tap = 0. (5) 

The corresponding diagrams are shown in Figs. 2 and 3, which have already been discussed 
in [12]. 

There has been some more work on the Black Box theorem: Hirsch et al. proved a su- 
persymmetric version of it [13, 14]. This extends the relation of the theorem to the lepton 
number violating scalar neutrino mass. We will not discuss that work further, as it is not of 
great importance for the discussion presented here. 

3 The Neutrino Mass generated by the Black Box Diagram 

Now we want to calculate the mass correction ("Butterfly") which is induced by the Black 
Box diagram shown in Fig. 4. We will first give a general parameterization of Qv(3(3 in 
terms of effective point-like operators, and then calculate the diagram for two examples of 
such operators in dimensional regularization. We will find that the Schechter-Valle theorem 
is not necessarily as useful as it might seem: The mass that is generated by the diagram 
in Fig. 4 is many orders of magnitude smaller than what is expected for the neutrino mass. 
Moreover, there exist operators generating but giving zero contribution to the Majorana 

neutrino mass via this particular diagram. Of course other diagrams, although more strongly 
suppressed, may give a non-zero contribution also for these operators. 

3.1 Parameterization of Neutrinoless Double Beta Decay 

The most general Lorentz invariant Lagrangian which contributes to 0^/3/3 via point-like 
operators has the following form [15]: 

C = -fm^^ {eiJJj + £2 J^V^^i + e^J^J^j + e^J^J^.f + e^Ji'Jj^) , (6) 

where Gp is the Fermi coupling constant, and rup is the proton mass. The hadronic currents 
required are given by 

J = n (1 ± 75) d, = (1 ± 75) d, J'^^ = u'- [7^ 7^ (1 ±l^)d, (7) 




Figure 4: The "Butterfly" diagram to be calculated in Sec. 3. 
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and the leptonic currents are 

j = e (1 ± 75) e^ f = (1 ± 75) . (8) 
In the remainder of this paper we will often use the abbreviations 

Pl = ^(1 - 75) and Pr = ^(1 + 75) (9) 

and write 

JL = e (1-75) e^ etc. (10) 

For all currents, different chirality structures are permitted (cf. Sec. 4). Note that in 
Eq. (6) we have already dropped the following terms which are Lorentz invariant, too: 

C = ^m-' {eeJ^rj^^u + ejJJ^-'j^, + eg J'^^j^) , (11) 

where the leptonic tensor current is given by 

r = e'-[j^,Y]{l±75)e''. (12) 

These terms had been included in the Lagrangian density in [15], but were neglected in the 
final analysis, as the authors worked in the s-wave approximation where these contributions 
vanish. In [16], it was pointed out that all operators proportional to 

t % 
nt^e", e- [7^, 7^] e^ and 675- [7^, 7^] e" (13) 

vanish identically because of the electron fields being Grassmann numbers. Therefore, the 
terms in Eq. (11) are not relevant for neutrinoless double beta decay. 

We may have to distinguish different chiralities of the currents involved, as in some cases 
the decay rate depends on them. Therefore, the authors of [15] used the parameters e^^^ 
with x,y,z = L/R to define the chiralities of the hadronic and leptonic currents in order of 
their appearance in Eq. (6). A suppressed chirality index indicates that it is not necessary to 
distinguish different chiralities. The limits obtained in [15] for the parameters are given in 
Tab. 1. We re-calculated the values of interest for us, using more recent results for the nuclear 
matrix elements {M.gt,n = 1-70 x 10~^ and M.f,n = 6.80 x 10~^) calculated in the Self- 
consistent Renormalized Quasiparticle Random Phase Approximation (SRQRPA) [17, 18]. 
Moreover, we have used qa = 1-25 and the half-life limit > 1.8 x 10^^ y obtained in the 
Heidelberg- Moscow experiment, as used in [15]. Note that, for 63, in addition to the updated 
values of the nuclear matrix elements, our cross-check of the bound from Ref. [15] resulted in 
a number that was smaller by a factor of 2 [19]. 

Let us mention that of course long range contributions to Oi^/3/3 — not proportional to 
the Majorana neutrino mass — may exist. A calculation of these diagrams, however, is not 
expected to add new insights. 
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k2| 


\^LLz\ \^RRz\ 
1' 1 


1 ^LRz 1 1 ^RLz 1 
1^3 1' F3 1 


|e4| 




from [15] 


3 X lO^'^ 


2 X 10^9 


4 X 10-*^ 


1 X 10"^ 


2 X 10"^ 


2 X 10"^ 


our calc. 


2.0 X 10-'^ 




1.5 X IQ-^ 









Table 1: Upper bounds for the coupling parameters in Eq. (6). They were evaluated "on 
axis," meaning that all other contributions were set to zero to extract the limits on one of 
the parameters. Different chiralities in the hadronic currents lead to different values only in 
the case of £3. 



3.2 Vertices and Propagators 

The Butterfly diagram we have to calculate (Fig. 4) is a non-standard one. At the effective 
vertex, lepton number is violated and two electrons are produced, which then leads to two 
outgoing neutrinos. We want to have a continuous fermion line in our diagram, so we have 
to rewrite some of the outgoing fields as incoming charge conjugate fields. The leptonic part 
of the standard electroweak vertex is given by 

-j= {c^PlvW^ + H.c.) . (14) 

We are interested in the Hermitian conjugate part, which is responsible for the annihilation 
of an incoming electron and the creation of an outgoing neutrino, and can be written as 

{e-i>'PLvW^)^ = V-f^'PLeW^ . (15) 

To rewrite it in terms of incoming charge conjugate fields, we have to transpose this expression. 
Using C7^'^C-^ = -7^, C'y^C-'^ = 75, and = -ip^C'^ we obtain 

{Vj'^PLef = ^i^Prv^ . (16) 

We additionally have to calculate the propagator of the charge conjugate electron fields: 



e%)?'=(x) . (17) 



Starting from the usual electron propagator 



with 



e{x)e{y) = iSpix - y) (18) 



5Hx-y)= / + r^. e-^(-^), (19) 



it is easy to see that 

e^{y)e^{x) = iSl{x-y). (20) 
If we now rewrite the expression in Eq. (17) with the help of 

e\y) = Ce^'iy) (21) 
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and the relation 

CSj^{x-y)C-' = SF{y-x), (22) 
which follows directly from the relations used to arrive at Eq. (16), we obtain 

e^(^{x) = Spiv - x) . (23) 

In the remainder of this section, we will work in momentum space. Therefore, a look at 
Eq. (19) reveals that the change of sign in the argument of Sp changes the sign of the 
momentum We thus have in momentum space 

' ZL — '6 + 777. 

e^{y)e'^ix)^i , ^+ . (24) 

3.3 Decay mediated by the Operator JlJlJl 
Our first explicit example will be the operator 

JlJlJl = SuPLcluPLdePLe" . (25) 

To calculate the diagram in Fig. 4, we have to write down the matrix element. Let us define 
the weak leptonic current 

jr=I77°Pz.e (26) 

and the weak hadronic current 

4 = dj^'PLU. (27) 
We then have to find all possible contractions in 

iflJLJddj'^W.^W^ffWaJ^W^li) , (28) 

where \i) and |/) denote initial and final state, respectively. It is easy to see that all contrac- 
tions lead to the same diagram given in Fig. 4. We find 



{f\i-f''PL^PLf^''-f''PRj ^dd-f^'PLUuPL^^^ . (29) 



Different contractions lead to 8 diagrams (there are 4 ways to contract the quark fields, and 
independently there are 2 ways to contract the W boson fields). Using the information from 
the last subsection and including all necessary factors, we can directly write down the matrix 
element of the process: 

f dhdhdq,dg,- -''^^ 



\k\-Ml){kl-Ml] 



j ' ' ---u. p^. „^ . ■■■^g ^p Pi Tr i ' " Pli " 7^ Pl 

{ki + qiY - ml qj-ml J V (^2 + 92)^ - ml - mj 



^yP)^ TTyi _^i.\2 2^ ' (30) 

[P — ^2) — TTlg [p + Kl)^ — ITlg 
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where we have used the short-hand notation dk = d^k / . The traces can be calculated 
in the standard manner, yielding 

V (^1 +qiY - ml qf-mj J [{ki + qi)^ - ml] {qj - m^) 

and 

\ik2+q2y-'ml qj-mj J [(^2 + 52)^ - m-2](g2 _ m^) 

We may simplify the last line of Eq. (30) by moving the Pls to the right. We thus project 
out the terms containing mg. The final expression for the self-energy of the neutrino at this 
order in perturbation theory is 

Y.(v) = — — ^ ' " ^ X, (33 

where the integral T is given by 

X = Xi X X2 , (34) 

with 



and 



Xi / dk.dq, ^ ^^^2 _ ^^^^^^^ ^ ^^^^ _ ^^^^^^ _ _ (35) 

^^2^9~2 _ ^^^2 _ ^2]p^ + ^2](^2 _ ^2)(^2 _ M^) ■ 

Integrals of this kind typically arise in two-loop calculations, which is not surprising: Our four- 
loop diagram factorizes because of the point-like nature of the corresponding operators, and 
it basically consists of two two-loop diagrams. These integrals may be calculated analytically 
in terms of Spence functions, for which extensive literature exists. Such calculations have, 
for example, been done in [20-22], in the framework of Higgs physics. For the moment, it 
suffices to know that the integrals are well behaved. To directly extract the dependence on 
masses and momenta, we will go another way: We can calculate Xi and X2 separately with the 
help of Feynman parameters and dimensional regularization. For example, let us calculate 
the integral 

Xi = / dq, , I[, (37) 

J {qf - m^) 

where 



I "^^^ [{p + k^Y - ml][{k, + q^Y - ml]{kl - M^) ' ^^^^ 



Using Feynman parameterization, we can rewrite 

-l-x 



:[= [ dx [ 

Jo Jo 



dy 



dki ^. (39) 

{x[{p + k,y - ml] + yPi + gi)2 -ml] + {l-x- y){kl - M^,)}^ 



The denominator can be rewritten as {[ki + {xp + yqi)]"^ — ^i}"^ , with 

A[ = {xp + yqif - xp^ -yql + {l-x- y)M^ . (40) 



8 



Here, we have neglected rrig and in comparison to M^. Performing the substitution 
fci — )• /ci — {xp + yqi), plugging the result into Eq. (37), and introducing a third Feynman 
parameter, we obtain 



167r2 
X / dqi 



1 rl-x 
dx 

^0 



dy I dz 



2 • 



{(1 - z){qf - ml) + z[{xp + - xp"^ - yqf + {1 - x - y)M^]} 

After completing the square and neglecting <^ Mw, we can write the denominator as 

2 



(41) 



{1- z + zy - zy) 



qi + 



xyz 



1 — z + zy^ — zy 



Ai{x, y,z) 



(42) 



where 

Ai(x,y,z) 



xyz 



{1 — z + zy"^ — zy) 



p 



(1 — z + zy"^ — zy) 



[zx'^p'^ — zxp^ + (1 — X — y)zM^ — m^] . (43) 



Plugging this result back into the expression for Xi in Eq. (37), performing one more linear 
substitution, dropping all terms linear in qi (as they vanish in the integration), and finally 
performing the integral over qi, we arrive at 



1 x 

dx 

^0 



dy 



X I dz 



xyz 



{1 — z + zy"^ — zy)^ \ € 



logAi-7 + log47r + C'(e) . (44) 



In this expression, e = 4 — d and 7 w 0.5772 is the Euler-Mascheroni constant. Calculating 
the integral I2 goes along the same lines. We finally obtain 



J: 



P 



(167r2)4 

X dx 
Jo Jo 

X da 
Jo Jo 



dy / dz 



xyz 



{1- z + zy^ - zy)^ \e 



- - log Ai - 7 + log 47r + ©(e) 



1 

db I dc- 



abc 



(1 - c + c62 - c6)3 Ve 
Ai is given in Eq. (43), and for A2 we have 



--logA2-7 + log47r + 0(e) . (45) 



A2(a, 6, c) 



abc 



{1- c + cb"^ - cby 



(1 — c + — cb 



[ca^p^ - cap^ + (1 - a - 6)cM^ - m^] . (46) 
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In the case of = 0, these expressions shnphfy to 



Ai(x,y,z) 



1 



((1 - X - y)zM^ - ml) 



(47) 



(1 — z + zy"^ — zy) 



and 



A2(a,6,c) 



1 



{{I - a - h)cMl, - ml) . 



(48) 



(1 - c + c62 - c6) 



The expression in Eq. (45) can be renormahzed via the usual "minimal subtraction" (MS) 
scheme. 

We are now able to write down an expression for the correction to the neutrino mass 
generated by the diagram in Fig. 4, assuming that there is no bare Majorana mass term in 
the Lagrangian. Plugging Eq. (45) into Eq. (33) we find 



where m^, is the physical neutrino mass. We have neglected the dependence on inside the 
logarithms. As the physical neutrino mass m^ is much smaller than M\y this does not affect 
the discussion. The physical mass m^ is the solution to the equation 



One solution of this equation is m^ = 0. The second solution m^ = a" , where a is the 
proportionality factor omitted in Eq. (49), is huge due to the smallness of a, and is unphysical 
as it does not correspond to any reasonable scale. Thus we obtain as physical correction to 
the zero neutrino mass: 



This means that the operator JlJlJl does not give a contribution via the Butterfly diagram 
shown in Fig. 4. We checked the calculations numerically with the TARCER program [23], 
which confirmed the zero result. 

One may now wonder if there exists an operator generating a non-zero contribution via the 
Black Box diagram at all. It does exist, and we want to give an example for such a diagram 
in the next subsection. The result, however, will be many orders of magnitude smaller than 
the expected light neutrino mass. 

Let us point out that the result of this subsection does not mean that all thinkable dia- 
grams involving the operator JlJlJl will give a zero contribution like the one we calculated. 
Any other diagram, however, will include more loops and will therefore be suppressed more 
strongly. Thus, the mass we expect to find will be even smaller than the one we find in the 
next section for a different operator. However, one should keep the result of this subsection 
in mind when arguing for the Majorana nature of the electron neutrino via the diagram in 
Fig. 4: Depending on the operator (i.e., on the underlying mechanism) of one possibly 

draws a diagram leading to a zero mass. Clearly, this does not support the hypothesis of a 
Majorana neutrino (although there might be other diagrams giving a non-zero contribution 
for the same operator). 




(49) 



(50) 



(5mi^ = . 



(51) 
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3.4 Decay mediated by the Operator J'^J^rJl 

In the previous subsection we have seen that the diagram calculated vanishes. However, there 
can of course be operators responsible for neutrinoless double beta decay which give non-zero 
contributions to the mass correction. In this subsection, we want to calculate the diagram 
assuming that the vertex is proportional to €3 from Eq. (6). Additionally, we choose the 
chirality structure in the following way: 

jRjRajL = SurPRdu^aPRdePLe" . (52) 

As in the former case, we have to find all contractions possible. So we can directly write down 
the matrix element of the diagram: 

f dhdhdg,dg, -''^^ 



w) 



V (^1 + 91)^ -ml qj-mj J \ (^2 + - ml " ql-ml J 

'^' \p-k,r-mf ' \p + k,Y-n.f '^^'^ • 
Calculating the traces in a similar manner as before, we obtain 

\\ki+qiY-ml^ ql-ml^ + qi? - mlU - m'^ ^^^^ 

and 

Tr i ^4^^^^^1.PJ-^,1^PL^ = -jr. ^^^^^^#1 T. ■ (55) 



{k2 + q2Y -ml qj-mj J [(/c2 + ^2)^ - ml] 

We can now rewrite the expression for the diagram as before: 



m 



di 



I28g-Gle,mlmlml^ ^ (56) 
nip 

where the integral I is given by 

T = Xi X X2 , (57) 

with 



and 



Xi / dk.dq, ^ _ ^^^^^^^^ ^ ^^^2 _ ^^^^^2 _ ^2)(^2 _ (58) 



^^2rf^~2 _ ^^^2 _ ^2]p^ + _ ^2](g2 _ ^2)(^2 _ ■ (59) 

The factor 4, which appears in Eq. (56) but not in Eq. (33), arises from a contraction 7'^7^ = 4. 
The numerators in Eqs. (58) and (59) do not contain loop momenta, which is the crucial 
difference to the case discussed in the previous subsection. 
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To calculate X\ and X2, we can combine the factors in the denominators containing k\ 
and ^2, respectively. Then, we may follow the calculations in [21], such that we arrive at the 
expression 

1 



Xi = X2 = ^^^2^^^4+, y ^((1 - zf^'^Mw, ma, m„; 0) 

with 

g{{l- z)^/^Mw,md,mu;0) = 

C2(z)e2 + ci(z)e + 7r^|^ + ^ [-1 + 27 + 2 ln(7r(l - z)M^)] + ^ + 



(60) 



1 7r2 
12 



+\ [-1 + 27 + 21n(7r(l - z)M^)]' - 1 + g{{l - z)'/^ Mw ,md,m^;0)^ 



(61) 



where the functions ci_2 are higher order corrections to the result from Ref. [21], which we 
have computed numerically. Note that these corrections can normally be neglected when 
computing ordinary two-loop diagrams. Since we, however, effectively need a product of two 
such diagrams, the higher order terms can appear in products with terms of order 0{e^^) 
or 0(e~^), thereby leading to non-negligible finite corrections. The function g appearing in 
Eq. (61) is defined by 



g{{l- zy/^Mw,md,mu;0) = 1 



U2 



Vl 



V2 



Ui 



— in a in ^= — 

2 ^ 



Sp I -r± ) + Sp ( -- ) + ^In^ - + ^In^ - + ^ln2 ^ - ^In^ ^ + ^ 



Ul 



Vl 



Vl 



where 



Mi,2 = 2 [l + b-a±VA' 

vi,2 = ^ {l-b + a±VA' 
A' = l-2{a + b) + {a-bf 



and 



(l-^)M^' 



(1 - z)M^ 



The Spence function is given by 



Sp(x) 



ln(l - y) 



y 



dy. 



, (62) 

(63) 

(64) 
(65) 

(66) 
(67) 



Integrating over z (the integration of 5 and ci^2 has been done numerically), and expanding 
in e, we finally obtain 



1 



(16^2)S 



C_2(M2,) C_i(M2,) 



+ Co(M^) + Ci(M^)e + C2{M^)e^ + O(e^) 



• (68) 
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Here, the coefficients are given by 

C-2{x) = 2, (69) 
C^i{x) = -3 + 27 + 21n(2;7r) -21n(47r^), (70) 

Cq{x) = 2.64493 + ^ [30 + 12(-3 + 7)7 + tt"^] + In(x^) [(-3 + 27) + ln(x7r)] 

+ 21n2(4^2) - [(-3 + 27) +21n(x7r)]ln(4^2), (71) 

Ci{x) = 4 - ^"^-fr ^ [30 + 12(-3 + 7)7 + + 12 ln(x7r) ((-3 + 27) + ln(x7r))l 
vr^ Iz 

- 2.64493 ln(47r2) +In2(47r2) (-3 + 27 + 2 ln(x7r)) -2 In3(47r2), (72) 

C^2(x) = 4-4l"(4^') 

+ ^ [30 + 12(-3 + 7)7 + 71"^ + 121n(x^) ((-3 + 27) + ln(x7r))] In2(47r2) 

+ 2.64493 In2(47r2) - \n^{AT:^) [(-3 + 27) + 2 ln(a;7r)] + 2 ln^(47r2), (73) 



and 



A= I dzci(2;) = 3.93180 X 10^ 

JO 

B= dz C2{z) = 2.28888 x 10^ . 



(74) 
(75) 



Then, 



f T 1 \ Cl,{M^^) , 2C_2(M^)C_i(M^) 
^ = ^^^^^=(16^1 ? + ^ 

2Co(M2,)C_2(M2,) + C\{Ml, 



+ 



+ 

2C_2(M^)Ci(M^) + 2C-i(M^)Co(M^) 
e 

+ C2(m2,) + 2C.i{M'^)Ci{Mlr) + 2C_2(m2,)C2(m2,) + 0(e)| 

^ (16^ t^o^^^/^'^ + 2C^,{M^^/^,')C,{M^^/^,') 

+2C.2{M^/f^^)C2{M^/fi^)] . (76) 

As we are interested only in an order of magnitude estimate for the neutrino mass cor- 
rection, we can renormahze the expression in Eq. (76) via the MS scheme (which means that 
we subtract the poles in e) and take the constant term as the leading contribution. Thus we 
find for the mass corrections 

^'^'^ ~ (167r2)4mp 

X [Cl{Ml,/n'') + 2C_i(M2,//i2)Ci(M2,//x2) + 2C-2{M^ I lJi^)C2{Ml, / ^?)\ . (77) 
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rriu 


= 1.5 to 3.3 MeV 




= 3.5 to 6.0 MeV 


rrie 


= (0.510998910 ± 0.000000013) MeV 


nip 


= (938.27203 ± 0.00008) MeV 


Gf 


= 1.16637(1) X lO-^GeV-2 


9 


= 0.652 


Mw 


= (80.399 ± 0.023) GeV 







Table 2: Masses and constants used in the calculation of Eq. (77) (taken from [24]). 



Using a renormalization scale = 100 MeV, as appropriate for a nuclear physics problem, we 
obtain 

6m^ = 9.4 X 10"^^ eV . (78) 

Here, we have used the values for masses and constants given in Tab. 2. It is obvious that 
this correction is far too small to be the main contribution to the neutrino masses we expect 
to have. In particular, such a mass cannot explain the neutrino oscillation data. Therefore, 
another mechanism must give the leading contribution of neutrino masses. However, if we 
know that a Maj orana neutrino mass is generated in a New Physics model anyway, we do not 
need the Butterfly diagram anymore. 

4 The Relation to Chiral Symmetry 

Finally, we want to consider the question in how far non-zero contributions to neutrino masses 
are related to the breaking of chiral symmetry. For example, in QCD it is well known that 
fermion mass terms necessarily break chiral symmetry, as they connect left-handed fermions 
with right-handed ones [25]. Hence, it is natural to ask the question whether the operators 
responsible for Oi^/3/3 violate chiral symmetry, as necessary for an operator that finally leads 
to a neutrino mass. 

Chiral symmetry acts on fermion fields as 

^p{x) e^"^5^(x) ~ (1 + ^075) ip{x) , (79) 

where the last step is valid for infinitesimal transformations. In this limit, it is easy to show 
that the usual right-handed Majorana mass term, which looks like 

^Majorana = -MrV'^-PrV' (80) 

in this notation, does not conserve chiral symmetry, as to be expected. The behavior of 
the different operators from Eqs. (7) and (8) is summarized in Tab. 3. Obviously, the only 
operators that are invariant under chiral symmetry are the ones with exactly one Lorentz 
index. Glancing at Eq. (6) we can easily see that, in order to transmit Oz^/3/3, the combinations 
of operators must be of the form "JJj," which is required by the simple fact that we need 
exactly two outgoing electrons. Such combinations of three currents can, however, never be 
built out of three operators with exactly one Lorentz index each, since a full contraction of 
the indices would never be possible. Accordingly, any operator that could be responsible for 
Oz^/3/3 must necessarily violate chiral symmetry. Turning this round, we have no operator in 
Eq. (6) that would be forbidden to yield a finite neutrino mass. 
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Operator Chiral Symmetry 



e (1 =F 75) not invariant 

67^^ (1 =F 75) e'^ invariant 

n (1 =p 75) d not invariant 

U7^ (1 =F 75) d invariant 

It|[7'^, 7^^] (1 =F 75) d not invariant 

Table 3: Behavior of the different operators under chiral symmetry. 



So why do not all operators lead to a non-zero correction? The simple answer to this 
question is that all of them will lead to some non-zero correction, but not necessarily at the 
order of the Butterfly diagram. The breaking of chiral symmetry is a necessary, but not a 
sufficient condition for the existence of the particular mass contribution under consideration. 
Accordingly, the Schechter-Valle theorem remains valid, with the only addendum that in some 
cases it might be necessary to draw a diagram that is more complicated than the one given 
in Fig. 1. 



5 Conclusions 

We have discussed in this paper the relation between the so-called Black Box operator (s), 
which are responsible for neutrinoless double beta decay, and neutrino mass operators. We 
emphasized that one has to be careful relating a possible observation of neutrinoless double 
beta decay fully to neutrino masses, while a very tiny Majorana contribution to the neutrino 
mass is guaranteed. The Black Box operator generates a Majorana mass contribution via the 
diagram in Fig. 4, which makes neutrinos Majorana particles, but we showed that the numer- 
ical contribution is many orders of magnitude too small to account for the observed neutrino 
masses. We have even identified one operator mediating 01/(3(3, but giving a zero contribution 
to the neutrino mass at leading (four loop) order in perturbation theory, demonstrating that 
the true Majorana contribution could even be further suppressed. We conclude therefore that 
lepton number violating operators beyond the Standard Model, which are not directly related 
to neutrino masses, guarantee radiatively non-zero Majorana neutrino mass terms, but the 
numerical values are extremely tiny such that other leading contributions to the physical 
neutrino mass must necessarily exist. If these mass terms are of Majorana nature, then they 
will dominate the Black Box diagram and there exists a direct relation to the observed rate of 
neutrinoless double beta decay. The Black Box induced Majorana mass is in this case only an 
extremely tiny correction found to be 5rn-^ = 0(10"^^ eV) which can for practical purposes 
be safely ignored. If, on the other hand, neutrinos are predominantly Dirac particles, then 
neutrinoless double beta decay is still possible due to other lepton number violating operators. 
The Black Box will then induce a tiny Majorana contribution, making neutrinos Majorana 
particles, but the numerical value is so tiny that one might consider this more or less as an 
academic statement. The important point is that neutrinoless double beta decay might exist 
as consequence of other lepton number violating physics beyond the Standard Model which 
is essentially unrelated to the values of the neutrino masses (an example for such a model can 
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be found in Ref. [26]). Let us rephrase this statement, to make it even clearer: Of course, if 
the neutrino is a Majorana particle and if we find that Oi//3/3 is predominantly mediated by 
light neutrino exchange, the experiments are able to test the neutrino mass scale. However, 
let us point out that we will not be able to find out if this is the case by experiments on 
neutrinoless double beta decay alone. This should be kept in mind when possible signals for 
neutrinoless double beta decay are discussed and translated into neutrino masses. 

Finally we would like to emphasize that our findings are logically fully consistent with 
the Schechter-Valle theorem, as the proof presented in Sec. 2 excludes that a Majorana mass 
term is forbidden by symmetry provided that neutrinoless double beta decay takes place. We 
have discussed how the well-known diagram may vanish at leading order, but at higher order 
there must exist a Majorana contribution. 
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